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^ ^^LZ^^li^^^y^n h^-^^ -oWent to *e -pie soMo^^^^^^^^ 

® o«H !^nmnlrtna thft TTV absorbanccs of the compound under reference conditions to die compound m a saturatea soiuuon, 
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TITLE OF THE INVENTION 
5 MEASUREMENT OF SOLUBlLITY-pH PROFILES 

CROSS REFERENCE TO RELATED APPLICATIONS 
This application claims priority to U.S. 
Provisional Application No. 60/178,616 filed January 28, 
10 2000, which is incorporated in its entirety herein. 

FIELD OF THE INVENTION 
The measurement of aqueous solubility in a high- 
throughput screening environment plays an important role 
15 in the selection of the most promising drug candidate 
molecules in pharmaceutical research and development. 
The invention described involves a method, a simple, 
robust, high-throughput screen, that is applicable for 
the determination of the equilibrium solubility of 
20 sparingly soluble compounds and that may be used an 
pharmaceutical, biotechnology, and related . industries . 
An analytical device has been designed to implement thxs 
solvibility measurement technique. 

BACKGROUND OF THE INVENTION 
The 1990 treatise by Grant and Higuchi [D.J.W. 
Grant, T. Higuchi, Solubility Behavior of Organic 
compounds, John Mley S Sons: New York, 1990] 
comprehensively covers the known art. Many protocols 
30 have been described in the literature [S. Venkatesh, J. 
Li, Y. Xu, R. Vishnuvajjala, B.D. Anderson, Pharm. Res. 
1996,13, 1453-1459; A. Avdeef, Pham. Pharmacol. Coimun. 
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«98 4. 165-178; A. J^vdeef, CM. Berger, C. Brown.U, 

.e.. .000, 17. .vdee* in B. .esta^ ■ . 

. r= — ^^^^^^^^^^^^ 

to Drug Permeation^ Absorptxorir Distribution, 

-jnni Ch. 22 (in press)] 
a.d Excretion, Onlv. Lausanne, 2001, < ^ 

for measuring ,olubillty-pH profiles, ■ usrng varrous 

detection systems. ,,,iled 
10 Classical approaches are based on the so-called 

satur.tion she.e-flas. method, and new rigorous methods 
,,„,Xly validated against it. However most 
..assical technigues are slow and cannot easUy be 
adapted for the high-throughput needs of modem drug 
15 discovery research, one focus of our invention. 

.t the early stages of research, drug candxdate 
compounds are stored as DMSO solutions, and .olubUxty 
measurements need to be performed on samples .n educed 
1„ OMSO, Often as 10-30 - solutions. It is welljno^ 
20 that even small quantities of MSO «5», .n water can. 
increase the solubility of molecules, and rt iS a 
Challenge to determine the true aqueous solub.Uty of 
compounds when DMSO is present, another focus of our 
invention. 

25 

Turbidimetric Ranking Assays ^ 

Turbidity •detection-based methods, popularized by 
Lipinski and others [A. Lipinski, F. Lonfcardo B.W 
Dominy, P.J- Feeney, Adv. Drug Deliv. Rev. 1997, 23 3- 
30 25; P.Q. Charmaine, M.B. Nicholas, A.K. Irwin, Eur. 
Phax.. Rev, 1998, 3(4); CD. Bevan, R.S. Lloyd, Anal. 
Chem. 2000, 72, 1781-1787], in part have met- some hxgh- 
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...oughput nee.s of d^, aUco,e.v — 

methods, aunou,. not — ^^^^^^'tr .,...ea 

attest to ran. molecules accord-, o 

i„hllitles. various inpiementations of tt. 

, are practiced at several phar.aceut.cal 

5 .ethod are P ^^^^^^ oetectron 

co.pan.es us.ng ^ ^^^^ „,p,elo.eters 

systems based on 36 veil labSystems: 

Franklin, MA, USA). The detector requires 

„ ........ -'»;t:s:ri t^»- — " 

-jr."- ■ ■rrr.i.; 

■ »^ uav The shortcomings of the 
customzed way. ^ ^^,„,i^iuty (in part due to 

methodology are P o repr ^^^^^^ 

" variabilrty of catt .i g ^^^.^....o, the 

sc.sti.es erratic -^^-^^ „alls 
sample suspension, and adhesion f 

, > the use of excessive amounts of DMSO, 
of vessels), the use ^^^^^ 

view-point standardisation of 

pharmaceutical industry, the lack or 



practice 

EFLC-Based Assays 



25 



30 



Based Assays 

. n ^«w.^=ir>if:.^ have taken the 
several pharmaceutical companies nave 

classicir saturation sha.e-flas. method and transferred 
it onto 96-weU plate technolo^ and a robotic li^id 

Iplsin. system, in an effort to --e «ut. 
^.alyses are performed "^p.op ate 

it is necessary to develop the app P 

^i,nric since in discovery, compounds 

chromatographic methods, since xn 

„.ay not be sufficiently characterized at the early 

^onpric fast gradient methods may be 
stages. However, generic tasi; gi: 
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. . .-n. the need for method development. In some 
eliminating the neea i.^.^^ted by a freeze- 

r.A^^ the DMSO is first eliminated Dy 
companies, tne uu buffers are added, 

drying procedure, before the aqueous buffers ar 
drying p ^ ^^^i can be 

which adds significantly to the assay ^^^^^^.^^ 

n Qamnles. Chromatographic 

problematic with volatile samples. 

Lection systems, although rich in J/^^J^ 
.erial, are inherently slow. Data handling is often 
rate-limiting step in the operations. 

Ideal solubility - PH relationships ,,„here 
Por ionizable molecules in dilute -^^^-- J"^;^;. 
salt precipitation may be ignored), when the log 
logarithm) Of the measured solubility is plot^e 

K< linear with one segment having a 
T^H the curve is bilinear, wj-uh 

, :;!.ope .na t.e o..« se^ent Ha..n, . u^t slop 

an acid such as diclorenac, 
example, for an acia suv- . . i araelv 

nH < 3 (where the acid is largely, 
is at a minimum for pH < J ^w"" 
L.a.,ea,. stays nearly constant .zero slope se^ent s 
,0 pH is increased, until the it exceeds the pK. (3.99), 
L: Which the solubility increases .y a decs e or 
every increased unit of pH, as more and .ore of the 
: : und converts fro. the uncharged state to the » - 
.olLle ne,atively-char,ed state (unit slope se,.ent^ 
« a Simple aqueous e^ilihriu. system, the pH where the 
. .L of the horizontal and diagonal segments 
:::::::: islgual to the p. of the somte. Onder s c. 
ideal circumstances, it is possihle to deter^ne the p^ 
a„d the intrinsic agueous solubility, 

vv,.. ornssina of the segments [Z-i. 
30 compound from the <=— ^ 1257-1260; «.H. 

Chowhan,, J. Piai»- „ ,,5, 

Stren,,- H.G.H. Tan, Int. a. PharH. 25. 135-145) . 
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s,eh pK. values ca. be directly measured by 
co«.er=lally-avaUable devices [CD. Bevan, 
D P Reynolds, Patent Cooperation Treaty, «0 99/X3328, 
18 March 1999; Sirius toalytlcal Instrument, Ltd., DKJ . 

complications »hloh May Ih.art tie Rellaile Measurement 

of Aqueous Solubility ^ 

certain surface-active compounds when dissolved rn 
„ater under conditions of saturation form self-. 

,0 associated aggregates or micelles, which can interfere 
with the determination of the true aqueous solubrlxty 
and the pK. of the compound IT.J. Roseman, S H. 
Valtows^y, ^. Pharm. Scl. 1973, «, 1680-1685,- C. Zhu, 
W.H. Strong, Int. J. Pharm. 130, 159-1681 . m^en 

15 the co.^.ounds are very sparingly soluble in water 
excipients are often added to enhance the rate of 

4-v,«> or-«>c!Pnce of the excipients 
dissolution. However, the presence 

also can interfere with the determination of the true 
aqueous solubility. If solubility measurements are done 
20 in the presence of simple surfactants ''^^.'''Z' 
Oh, J.R. Crison, G.L. Amidon, J. Pham. Sea. 2000, 89, 
268-274], bile salts [S.D. Mithani, V. Bakatselou C . 
TenHoor, J.B. Dressman, Phana. Res. 1996, 13, 163-167], 
cyclodextrins [P. Li, S.E. Tabibi, S.H. Yalkowsky, J. 
25 Pharn. Sci. 1998, 87, 1535-1537], or ion-pair formxng 
counterions [J.D. Meyer, M.C. Manning, Phar.. Res, 1998, 
15 188-1931, extensive considerations need to be 
applied in attempting to extract the true aqueous 
solubility from the data, third focus of our invention. 
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BRIEF SUMMARY OF THE INVENTION 
- - - This invention relates to a method for the 
determination of solubility of a compound and an 
analytical device for carrying out said method. The 
basic method involves determining solubility of a 
compound by measuring the UV spectrum of a reference 
solution of the compound, under conditions avoiding or 
suppressing precipitation, and comparing it to the UV 
spectrum of a saturated sample solution of the compound, 
variations of the basic method include: (a) making 
reference solutions either by dilution of the sample 
solution to the point where precipitation is avoided, or 
making reference solutions by adding a water-miscible 
cosolvent to the sample solution so that precipitation • 
is suppressed, and comparing the UV absorbances of the 
compound under reference conditions to the compound in a 
saturated solution, (b) determining the true aqueous 
' solubility from the effect on the pK. that results from 
dissolving the compound in an aqueous solution 
containing some DMSO (typically 0.1-5% v/v) , and (c) 
correcting concentrations determined from the UV 
absorbance values for impurities and other factors that 
rnight affect the shape of the sample absorbance curve 
taken of the saturated solution. 

solubility determinations are illustrated with 15 
sparingly-soluble generic, mostly-ionizable drugs (Fig. 
1), covering about- three orders of magnitude xn . 
solubilities, with the lowest value being about 0.1 
j.g/mL (terfenadine at high pH) . The ionizable bases 
30 include amiloride, amitriptyline, chlorpromazine, 
miconazole, nortriptyline, phenazopyridine, propranolol, 
and terfenadine. The ionizable acid ^ set consists of 
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diclofenac, furosemide, indomethacin, 2-naphthoic acid, 
and probenecid. Piroxicam was picked as an example of 
an ampholyte, and griseofulvin as an example of a non- 
ionizable molecule. The results of the solubility 
measurements were compared to literature equilibrium 
solubility-pH values from reliable sources, and to ' 
values determined by the pSOL Model 3 (pION) acid-base 
titrator [A. Avdeef, CM. Berger, C. Brownell, PJiarm. 
J^es. 2000, 17, 85-89]. 
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BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 
Figure 1 shows the structure of compounds evaluated 
using the invention and for which data appears herein. 

Figure 2 represents a schematic of an analytical 
device utilized to carry out the methods described in 
this patent. 

Figure 3 is a flow chart showing the aqueous 
dilution method: (a) general depiction, (b) specific 

implementation . 

Figure 4 is a flow chart showing the cosolvent 
method: (a) general depiction, (b) specific 
implementation . 

Figure 5 illustrates the solubility - pH curve for 
a hypothetical ionizable acid with a pKa of 4 and three 
25 different intrinsic solubility values: (a) 5 jig/mL, (b) 
50 jig/mL, and (c) 250 ^g/mL. In each of the cases, the 
upper horizontal line represents the maximum possible 
sample solution concentration, and the lower horizontal 
line represents the maximum possible reference solution 
concentration, of. Fig. 6 below. Intersections of the 
lines with the curve form the working range, of the 
aqueous dilution method for ionizable compounds. 
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Figure 6 shows the expected ratios of the areas - 
under-the-curve, R (equation 5 below), corresponding to 
the examples in Fig. 5. Solubility values are derived 
under the conditions 0.05 .< R < 1 (equation 6). Note 
5 that 0.05 refers to the 20-fold dilution ratio in the 
aqueous dilution method. 

Figure 7 shows curves generated from data, using 
equation 6, for solubility as a function of pH for 
several of the compounds studied by the aqueous dilution 
10 method. The horizontal lines represent the maximum 
possible sample concentration, Cs. The sample solution' 
is not saturated when the discrete points coincide with 
the horizontal line, indicating that the solubility 
cannot be determined for those points. 
15 Figure 8 shows for several of the studied molecules 

the logarithm of the apparent solubility versus pH 
curves in solutions containing 0.5% (v/v) DMSO (solid 
line), and the corresponding calculated aqueous 
solubility versus pH, using equations from Table 1 below 
20 (dashed line). The circles represent measured points; 
the solid line represents the best-flit curve of the 
measured points, in accordance with the expected shape 
of a log S vs. pH curve; the dashed line represents the 
log S vs. pH curve after correction for the anomalies ■ 
25 caused by mSO, aggregation, and the like. 

Figure 9 shows UV spectra of piroxicam taken at 
several values of pH, typical of spectra taken in the 
aqueous dilution method. Note that the reference 
(right) absorbances are 20 times lower in optical 
/30 densities, compared to the highest absorbances of the 
sample spectra (left) . 



- 8 - 



wo 01/55698 



PCT/USOl/02377 



Figure 10 shows UV spectra of (a) miconazole and 
(b) phenazopyridine, taken at jseveral value's of pH, . 
typical of spectra taken in the cosolvent method. Note 
that the reference (right) absorbances are slightly 

5 higher than the highest sample absorbances (left) , in 
contrast to the situation in the aqueous dilution method 
(Fig. 9). In principle, this indicates that the 
cosolvent method is more sensitive than the aqueous 
dilution method. 

10 Figure 11 shows sample and reference spectra, where 

there is a shape difference in the region of 250-350 nm. 
The special weighting scheme, equation 27 below, down- 
weights the data in the region of 285 nm by a factor of 
169. Multiplying the reference curve by 10 and 

15 superimposing it on the sample curve produced the dashed 
curve . 

Figure 12 shows the alkalimetric pH titration curve 
of the aqueous universal buffer solution. 

Figure 13 shows the alkalimetric pH titration curve 
20 of the cosol vent-water universal buffer solution. 

DETAILED DESCRIPTION OF THE INVENTION 
We have adapted the classical shake-flask method to 
the 96-^311 microtitre plate format, but our detection 
25 system is not based on the traditional, one-compound-at- 
a-time HPLC method. The invention is a parallel 
detection method, based on the use of a 96-well 
microtitre pl^te UV/vislble (hence called 'UV') 
spectrophotometer. Our invention" can measure 96 
30 solubilities at a time, significantly increasing 
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throughput. Compound-dependent method optimization is 
not required in the direct DV method. 

(1) The essence of our invention is that of making 
concentration standards under conditions avoiding or 
5 suppressing precipitation for IN spectrophotometers (in 
two different implementations: (a) aqueous dilution and 
(h) cosol vent; , using molecules, especially ionizable • 
molecules, whose UV spectroscopic properties are unknown 
at the start of the assay. (Simply weighing the 
10 compounds and preparing solutions does not work in 
general, because the compounds precipitate under 
conditions necessary to characterize their 

solubilities. ) 

(2) Also, a novel and simple computational method 

15 is developed to extract the aqueous intrinsic 
solubilities of drug molecules from data altered by 
DMSO-drug binding, or drug-drug aggregation reactions, 
or drug reactions with non-ionizable excipients, using 
the difference between the apparent pK, measured by the 

20 analytical device and the true pKa measured by . 
analytical devices available commercially (e.g., from 
Sirius Analytical Instruments Ltd., Forest Row, E. 
Sussex, UK) . 

(3) Also, an improved general method for 
25 determining ' concentration by OV spectrophotometry is 
derived, based on using a novel peak shape algorithm for 
adjusting weights in a whole-spectrum weighted 
regression analysis, matching spectra of reference 
solutions (of known . concentration, under, conditions 
30 avoiding or suppfessing precipitation) 1:c solutions 
containing analyte of unknown concentration (due to 
precipitation) . 

- - 11 - - 
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DEFINITIONS OF TERMS 
If the molecule is ionizable, its charge depends on 
the pH of the solution in which it is dissolved. The pH 
5 at which the concentration of the charged form is equal 
to the concentration of the uncharged form is called the 
pKa of the molecule. For example, for the ionizable 
acid (or ionizable base), represented by HA (or B) , 
chemists usually denote the equilibrium reaction by the 
10 equation 

HA ^ A- + ( or BH^ ^ B + hM (D " 

and the corresponding equilibrium ionization constant is 
15 defined as 

Ka - [A-][H*] / IHA] (or = [B] [H*] / [BH"^] ) (2) 

By definition, pKa = -log Ka. VJhen a solute molecule, 
20 HA (or B) , is in equilibrium with its precipitated form, 
HA(s) (or B(s) ), chemists usually denote the process 
by the equilibrium expression 

HA(s) ^ H (or B(s) ^ B) (3) 

25 

and the corresponding equilibrium constant is defined as 

So = [HA]/[HA(S)] = [HA] (or So = [Bl/[B(s)] = [B]) (4) 

/ 

So By convention, [HA(s)] = [B(s)] = 1. Equations (3) 
represent the precipitation equilibria^ of the uncharged 
species, and are characterized by^ the intrinsic 

- 11 - 
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solubility equilibrium constant, So- • (The zero 
subscript denotes -the zero charge of the precipitating 
species . ) 

The values of the equilibrium constants in • 
5 equations (2) and (4) depend on the solvent system used 
(water and aqueous buffers in our case), the inert salt 
used to adjust the ionic strength of the solution (e.g., 
NaCl or KCl) , temperature {22-25°C in our case), 
pressure (ambient in our case), and may be affected by 
10 organic components (other than the sample) in the 
aqueous solution, such as DMSO, propanol, surfactants, 
bile acids, lipophilic counter ions, cyclodextrins, and 
the like. For example, if the solvent system were 
changed from water to 50% v/v methanol-water, then the 
15 pKaS of acids would be higher and of bases would be 
lower, compared to values determined in pure water [A. • 
Avdeef, J.E.A. Comer, S.J. Thomson, Anal. Chem. 1993, . 
65, 42-49] . 

We use the term solubility to mean the 
20 concentration of a solute (in units of moles per liter, 
M, or micrograms per milliliter, jig/mL) in a saturated 
pH-buffered aqueous solution. Under pH conditions where 
the sample molecule is essentially uncharged, the 
measured solubility is equal to the intrinsic solubility 
25 of the molecule. Solubility may depend on pH, but the 
intrinsic solubility, being a thermodynamic equilibrium 
constant, does not. 

We will implicitly refer to pKaS or log So or log S 
of compounds as values determined in water or an aqueous 
30 buffer. ^' When we use the term "apparent" pK^, or 
"apparent" log So, we do so to emphasize that ^the 
aqueous value may not be evident when solutions contain 
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components such as DMSO, propanol, surfactants, bile 
acids, lipophilic counterions, cyclodextrins , and the 
like, that may affect the ionization properties of the 
compound. 

5 The term "stock" solution will refer to a solution 

made from a precisely known quantity of pure compound 
dissolved in a known volume of solvent, usually pure 
DMSO, usually at a 10-100 mM concentration level. 
Compounds are expected to be fully dissolved in the 
10 stock solution. 

We use the terra "reference" solution to mean a 
solution of the sample compound where the concentration 
of said compound is known. The reference solution can 
be prepared by precisely weighing a quantity of the pure 
15 compound and dissolving it 'in a precise volume of 
solvent (usually water or aqueous buffer) , or by adding 
a precise volume of a stock solution of said compound to 
a precise volume of solvent, under conditions avoiding, 
or suppressing precipitation. Two ways are developed to 
20 suppress or avoid precipitation: (a) simply make a 
dilute aqueous solution, below the limit of solubility, 
or (b) add water-miscible cosolvent to an otherwise 
saturated aqueous solution. 

The term "saturated" solution refers to a solution 
25 at equilibrium, containing enough sample compound, so 
that' a portion of the compound is precipitated out of 
the solution and a portion remains dissolved, at a 
concentration equal to the solubility. In other words, 
a saturated solution of a compound contains an excess of 
30 undissolved solid form of the compound in contact with 
its solution. 
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The term "supersaturated" solution refers to a 
solution containing no solid and more dissolved compound 
than expected from its equilibrium soltability. Usually, 
such solutions are unstable, and given enough time, will 

5 precipitate . 

By the term "filtration", we broadly mean the 
separation of the precipitated solid from a saturated 
solution, either by passing said solution through a 
filter, by sedimentation, by centrifugation, or by any 
10 other separation means. 

Saturated solutions are usually "concentrated," and 
reference solutions are usually "dilute," unless 
cosolvent is used. 

To cause the formation of a saturated solution, a 

15 small aliquot of the stock solution (typically 1 - 50 p-L) 
or a weighed quantity of compound is added to an aqueous 
buffer solution (typically 1000 - 2000 pi), which is 
called the "system solution" here. 

By the term "excipient", we mean any additive in a 

20 solution that may affect the solubility of the solute. 
More specifically, we mean a non-ionizable additive in 
the pH interval where the sample compound is ionizable. 
Examples of such additives, as we use the term, are 
DMSO, surfactants, bile salts, phospholipids, 

25 cyclodextrins, ion-pair forming count erions, polymer s^. 
cosoi vents, and combinations thereof. 

COMPONENTS OF THE ANALYTICAL DEVICE OF THE INVENTION 

/ 

30 Hardware 

The analytical device, shown in a block diagram in 
Fig. 2, consists of a robotic liquid handling system 

- 14 - 
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(1) , a raicrotitre plate scanning UV spectrophotometer 

(2) , a pH titrator device {3), a microtitre plate vacuum 
filtration manifold (4), a microtitre plate washer (5), 
a microtitre plate orbital shaker {6), four precision 
syringe dispensers ( 7) , four dispenser arms positioned . 
by the robot anywhere on the worktable of the liquid 
handling system (fl), a wash station and waste trough 
(9), a small-tips pipet rack holder (10), a large-tips 
pipet rack holder (11), a used-tip collector (12), a 
stock sample microtitre plate {13), a 20-fold diluted . 
stock sample microtitre plate (14), a deep-well 
microtitre plate for saturated solutions (15), a deep- 
well microtitre plate for reference solutions in the 
aqueous dilution method {16) , a plastic UV microtitre 
plate (17) , four test tubes filled with rinse solution 
(IS), a test tube for cosolvent (15), a test tube for 
the DMSO {20), a test tube for 0.5 M KOH titrant {21), a 
titration vessel, with a magnetic stir bar and a 
magnetic stir motor underneath (22) , a test tube to 
store the pH electrode . {23), and an electrode wash 
station (24) . 

The robotic liquid handling system (1) , is equipped 
with four high-precision 0.5 mL syringes equipped with 
two-way valves (7), and four dispenser arms, with 
25 conductive-plastic disposable tips (8), a liquid-level 
sensing system, two racks for tips (10,11), six 
microtitre plate holders (13-17), dram-size vial and 
test tube holders (18-21), a rinse and waste station 
(9), and approximately 100x50 cm workspace (e.g., 
Genesis 100/4 System from Teean, --Research Triangle Park, 
NC, USA ) . 



20 



30 



- 15 - 



wo 01/55698 



PCTAJSOl/02377 



A 96-well microtitre plate scanning 
spectrophotometer is a part of the analytical device 
(2), with wavelength range at least 200-850 nm, with <5 
nm resolution (e.g., SpectraMAX 190 from Molecular 

5 Devices, Sunnyvale, CA, USA) . 

Although the 96-well microtitre plate format is 
used preferentially, the method and the analytical 
device can be constructed with other plate formats, for 
example, the 384-well microtitre plate format. 

10 Although our preferred detector system is a . 

scanning microtitre plate UV spectrophotometer, diode- 
array microtitre plate UV spectrophotometers, and flow- 
through UV detector systems, scanning or diode array, 

may also be used. 

15 Although a UV detector is the preferred embodiment 

of the invention, the detector may be any suitable 
spectrophotometric analytical detector, such as 
ultraviolet or visible spectrophotometer, a fluorimeter, 
a colorimeter, a light-scattering device, polarimeter, 

20 optical rotation or circular dichroism detector. 

A pH titrator (3) is used, having a vessel to hold 
the solution being titrated (22), equipped with a pH 
meter capable of precisely reading pH from 1.5 - 12.5, 
having a dispenser able to add 0.5 M KOH titrant in 

25 precise small .amounts, such as 1 ^iL (7), and able to 
stir the solution during titrant additions (e.g., pSOL 
Model 3 from pION, Woburn, MA, USA) . 

A vacuum filtration manifold (e.g., from Whatman) 
for microtitre plates (4) is used, with a source of 

30 vacuum. 

A microtitre plate washer (5) is used, with a 75% 
v/v methanol wash solution (e.g., from Tecan) . 
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A microtitre plate orbital shaker is used {6) 
(e.g., from Lab-Line- Instruments, Inc.). 

Also used are commercially-sourced .96-well 
polyethylene/polypropylene microtitre plates, in 0.4 and 
5 2.2 itiL well capacities (numerous sources); coramercially- 
source plastic UV plates (e.g., from Coming-Costar or 
Greiner) . 

In the basic protocols described below, 8 compounds 
are sampled in each microtitre plate, at 12 different 
10 pH-buffer solutions. This reflects the 8-row x 12- 
column layout of the microtitre plates. The procedures 
can be easily scaled up three-fold to 24 samples per day 
on the present robotic system, still preserving the 12- 
pH profile, or to 384 samples per day at one pH. 

15 

.Reagents 

A universal buffer system solution is used (25), 
designed to have a linear response to alkali titrant 
additions in the pH range 3-10, with nonlinearity in pH 
20 of about ± 0.05 pH units. The solution possesses buffer 
components with low-UV absorption (each component with 
OD < 0.05 at 220 nm for a 4.5 mm path length). The 
solution possesses buffer components with low tendency 
to interact with sample species, and specifically 
25 excludes phosphate and citrate. A 2-L solution of the 
buffer, the pH of which is initially near 3, has the 
capacity of about 100 mL of 0.5 M KOH when pH is raised 
to 10. The ionic strength of the universal buffer 
solution is« about 0.01 M. Such a universal buffer 
30 solution i^s specifically designed for solubility 
measurements and is commercially available from pION. 



M - - 
I 
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A Standardized 0.5 M KOH solution is used (21), 
containing < 0.5% M/M carbonate 1 available from a number 
of commercial sources) . 

Organic solvents are used: spectroscopic grade (or 
5 better) DMSO {20) and 1-propanol (19) (numerous 
commercial sources) . 

Software 

A computer program is used, which controls the 
10 actions of the robotic fluidic delivery system, prompts 
the operator to perform certain tasks, controls the 
actions of the spectrophotometer, and processes the 
spectral data to determine aqueous solubilities as a 
function of pH and the intrinsic aqueous solubility, 
15 displays the data graphically and in report forms, and 
which transfers the results of the analyses to a 
Microsoft® Excel spreadsheet. It is commercially 
available from pION. 

The software is loaded on a Windows NT® computer, 
20 which connects to the robotic workstation (1), the UV 
spectrophotometer (2) , . and the pH titration device (3) , 
by RS232 serial cables, and controls the actions of the 
analytical device and stores the spectra collected for 
further processing. 



25 



30 



IMPLEMENTATION BASED ON AQUEOUS DILUTION METHOD 

The diagram in Fig. 3a shows the general concept of 
the aqueous dilution method of preparing reference 
solutions for UV measurements. The "sample" track on 
the -left side of Fig. 3a depicts a series of steps 
leading up to the taking of the sample spectrum. A 
known quantity of sample {27) is added to a known volume 
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of the universal buffer solution of known pH (28), the 
amount of said sample being sufficient to cause 
precipitation to occur in the formed saturated solution 
(25). After waiting a period of time to allow the 
5 saturated solution to reach a desired steady state, the 
solution is filtered (30) to remove the solid, to obtain 
a clear solution (31), whose spectrum is then taken by 
the UV spectrophotometer (32). The mathematical 
treatment of the spectral data (described below) yields 
10 the area-under-the-curve of the filtered sample solution 
(31), AUCs. 

The "reference" track on the right side of Fig. 3a 
depicts a series of steps leading up to the taking of 
the reference spectrum. A known quantity of sample (33) 
15 is dissolved in a known volume of the universal buffer 
solution of known pH (34), the amount of said sample . 
being X-times less than 27, in order to avoid 
precipitation in the formed solution (35) . The spectrum 
is immediately taken by the UV spectrophotometer (35), 
20 to take advantage of the possibility that solution 35 
may be "supersaturated" (that is, solid should have 
precipitated, but because not enough time was allowed 
for the solid to precipitate, the solution was 
temporarily clear and free of solid) . The mathematical 
25 treatment of the spectral data yields the area-under- • 
the-^urve of the reference sample solution (35), AUCr. 
When the ratio of the two AUC values, 

R = AUCr / AUCs 

30 / 

is greater than 1, supersaturation is indicated, and 

solubility can be determined by equation (6) below. 



- 19 - 



wo 01/55698 



PCT/USOl/02377 



10 



When R = 1, the sample and the reference solutions 
are both saturated, and solubility cannot be determined, 
unless visual inspection of solution 35 indicates that 
no solid is present, which could happen in a 
supersaturated solution. In this case, solubility can 
be determined by equation (6) below.. 

TO overcome the situation of R - 1 with both 
solutions being saturated, a higher value of X needs to 
be tried, and the procedure repeated until R < 1. 

■ If solubility cannot be determined because R - 1, 
one can conclude that the solutility is less than or 
equal to the calculated concentration of 3S (the 
calculation ignoring precipitation) . 

At the other end of the" scale of ratios, when R - 
1/X, the compound is fully dissolved in both the pre- . 
filtered sample solution (29) and the reference solution 
(35) , and solubility cannot be determined. (Ideally, R 
cannot be less than 1/X, unless Beer-s law is violated.) 
However, one can conclude that solubility is greater 
than or egual to the calculated concentration of the 
solution 35 ^concentration C.) . To overcome thxs 
particular limitation, more sample 27 may be used. 

When the condition 1/X < R < 1 is «iet, then the 
25 solubility of the sample compound is 

S - Cr / R 



15 



20 



30 



where Cr is the calculated concentration of 35. 

in example 1, we implement the aqueous dilution 
method with the value of X. = 20. 
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IMPLEMENTATION BASED ON COSOLVENT METHOD 

The diagram in Fig. 4a shows the general concept of 
the cosolvent method of preparing reference solutions 
for UV measurements. The "sample" track on the left 
side of Fig. 4a depicts a series of steps leading up to 
the talcing of the sample -spectrum in a water-cosolvent 
solution. A known quantity of sample (37) is added xn a 
known volume of the universal buffer solution of Icnown . 
pH (38), the amount of said ss^le being suffxcxent to 
, cause precipitation to occur in the formed saturated 
solution (39). After waiting a period of ti^e to allow 
the saturated solution to reach a desired steady state 
the solution is filtered (40) to remove the solid, and 
obtain a clear solution (41) . 
5 A volume Y of a water-miscible cosolvent 42 xs 

added to a volume Z of solution 41 to produce a new : 
solution 43, in which the compound is diluted by • 
Z/(y+Z), from its concentration in solution 41. The 
spectrum of solution 43 is then in^ediately taken by the 
10 UV spectrophotometer (44), to minimize evaporation of 
the water-cosolvent solution. The mathematical 

treatment of the spectral data yields the area-under- 
the-curve of the filtered cosolvent sample solutxon 
(43), AUCs~^ 

25 suitable cosolvents may be selected from a nu^er 

of - water-miscible organic solvents, such as 
acetonitrile, methanol, ethanol, iso-propanol, 1- 
propanol, ethylene glycol, propylene glycol, 
polyethylene glycol 400 (PEG-400) , 1,4-dioxane, 
30 dimethylformamide, acetone, tetrahydrofuran; dimethyl- 
sulfoxide (DMSO), and mixtures thereof. The ones wxth 
the lowest vapor pressure, the greatest capability xn 
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dissolving a solute (i.e., highest solubilizing power), 
and the lowest UV absorption are preferred. 

The "reference" track on the right side of Fig. 4a 
depicts a series of steps leading up to the taking of 
5 the reference spectrum in a water-cosolvent solution. A 
known quantity of sample .(45) is added to a known volume 
of the universal buffer solution of known pH {46) , the 
amount of said sample being comparable to that of 37, no 
effort being made in this step to suppress precipitation 
10 in the formed solution (47) . 

A volume Y of a water-miscible cosolvent 42 is 
added to a volume Z of solution 47 to produce a new 
solution 48, in which the compound is diluted by. 
Z/(Y+Z), from its concentration in solution 47. 
15 It is necessary to ensure that the cosolvent 

completely dissolves the compound in solution 48. If 
this is not the case, different values of Y and/or Z 
must be chosen and the whole procedure repeated to 
ensure that solution 48 is solid-free. 
20 The spectrum of solution 48 is then immediately 

■ taken by the UV spectrophotometer (49) , to minimize 
evaporation of the water-cosolvent solution. The 
mathematical treatment of the spectral data yields the 
area-under-the-curve of the filtered cosolvent sample 

COS 

25 solution (48) , AUCr . 

,.- Define the ratio of the two AUG values, 

rCos = aUCr°°^/ AUCs"^ 
30 The solubility of the sample compound is ' 

S = (1 + Y/Z) Cr"= / 
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Where C^"^ is the calculated concentration of the 
compound in solution 48. 

in example 2, we implement the cosolvent method " 
5 with the value of Z/(Z+Y) = 0.5, using 1-propanol. 

METHOD USING pK^SHIFTS FOR DETEmiNING TRUE AQUEOUS 
SOLUBILITY OF lOmZRBLE COMPOUNDS IN THE PRESENCE OF 
NON-IONIZABLE CHEMICALS THAT DISTORT THE SOLUBILITY-pH 
10 PROFILES 

DMSO Binding to the Uncharged Form of a Compound 

Although the results depicted in Fig. 7 (plots of 
experimental S vs. pH) were very satisfactory at first 
15 sight, we soon discovered unexpected anomalies, which 
are not immediately apparent from the figures. The • 
aberrations became evident when log S were plotted vs. 
pH (Fig. 8) . 

we found that the log S vs. pH curves (cf., 
20 Background section) were altered in the Presence of 0.5% 
v/v DMSO, in that the apparent pKa values, pKa^^ derived 
from log S vs. pH, were different .from the true pKe 
values by about one log unit. The pK^ values were 
generally higher than the true pKaS for ionizable acids 
25 ("positive" shift), and lo^er than the true pKaS. for 
ionizable bases ("negative" shift) . 

■ we believe that this is caused by DMSO binding to 
the drugs. Although we did not anticipate it at first, 
we actually had collected the data necessary to 
30 calculate the extent of dr^^g - DMSO binding, and also to 
- correct for the effect, to obtain true aqueous 
solubility -pH profiles. For example, equation (9) 
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represents the binding of the uncharged species by n 
molecules of DMSO (e.g., in 0.5% v/v) . 

HA + n DMSO ^ HA{DMSO)n • (9) 

5 

If the effect of the DMSO is present, and it is ignored 
in the analysis, then the determined solubility is not 
the value one would measure in the absence of DMSO. 

We discovered that if one determines the pKa*" (in 

10 the presence of the DMSO effect) in the solubility ' 
assay, and subtracts from it the true pKa (determined 
using commercial instruments such as those from Sirius, 
UK), the difference, when subtracted from the logarithm 
of the apparent (DMSO-distorted) solubility. So , 

15 yields the true aqueous solubility constant: 

log So - log So*^* - {pKa*** - PKa) r 

for ionizable acids (10) 

20 DMSO makes the compound appear more soluble, and the 
true aqueous solubility can be determined from the 
apparent solubility by subtracting the pKa difference, 
as done in equation (10) . 

Furthermore, we discovered that if one is 

■25 considering the case of an ionizable base, such as 
chlorpromazine, the shift of the pKa is to values loiter 
than the true pKa. The corrective treatment is 
indicated by equation (11) 

30 , log So"= log So*'^ + (pKa^^^ -^^pKa) , 

for ionizable bases (H) 
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For an amphoteric molecule (which has both acid and 
-base functionality) with two pK^s, either equation (10) ' 
or equation (11) may be used, depending on which of the 
two pKaS is selected. 

Fig. 8 illustrates the apparent solubility -pH 
curve (solid line) and the true aqueous solubility-pH 
curve (dashed line), correcting for the effect of DMSO 
for several of the molecules considered. 



/ 
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Uncharged Forms of Compound-Compound Aggregation 

Shifts in the pKa can also be expected if water- 
soluble aggregates form from the uncharged monomers. 
This may be expected with surface active molecules or 
5 molecules such as piroxicam [J. Jinno, D.-M. Oh, J.R. 
Crison, G.L. Amidon, J- Pharm. Sci. 2000, 89, 268-274]. 
Consider the case where no DMSO is present, but 
aggregates form, of the sort 

10 m HA 5± {HA)ni (12) 

The assumption is that the aggregates are water soluble, 
that they effectively make the compound appear more 
soluble, and if ignored, they will lead to the erroneous 
15 assessment of solubility. We discovered that equations 
(10) or (11)/ for ionizable acids or ionizable bases, 
respectively, applies also to the case of aggregation. 

Compound-Compound Aggregation of Charged Ionizable Bases 
20 Consider the case of an ionizable base, where the 

protonated, positively charged form self-associates to- 
form aggregates, but the uncharged form does not. We 
may have this case with phenazopyridine (Fig- 8). 
Phenazopyridine is a base which consistently shows 
25 positive shifts in its apparent pKa/ the opposite of 
what's expected of uncharged-compound DMSO/ aggregation 
effects. A rationalization of this effect can be based 
on the formation of partially protonated aggregates 
(perhaps micelles) . For example, one can hypothesize 
30 that one of the species is (BH'**)n. 

n BH'' 5± (BH*)n (13) 
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We discovered that for such a case, the obs.erved 
solubility-pH curve is shifted horizontally (not 
vertically, as with uncharged-compound DMSO/aggregation 
5 effects), and that the apparent intrinsic solubility is 
not affected by the phenomenon. 

Compound-Compound Aggregation of Charge lonizable Acids 
Consider the case of an ionizable acid, where the 

10 de-protonated, negatively-charged form, is strongly 
self -associated, but the uncharged form is not. We may 
have this case with prostaglandin F2a [T.J. Roseman, 
S.H. Yalkowsky, J. Pharm. Sci. 1973, 62, 1680-1685]. 
Prostaglandin is an ionizable acid with a reported 

15 negative shift in its apparent pKa, the opposite of 
what's expected of uncharged-compound DMSO/aggregation . 
effects. A rationalization of this effect can be based 
on the formation of partially deprotonated aggregated 
complexes (perhaps micelles) . One can hypothesize a 

20 reaction similar to equation (12)', but involving 
negatively-charged ions. The reported solubility-pH 
curve is shifted horizontally (not vertically, as with 
uncharged-compound DMSO/aggregation effects), and that 
the apparent intrinsic solubility is not affected by the 

25 phenomenon. 



/ 
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lonizable Compound Binding hy Non-ionizable Excipients 

We have concluded that a number of phenomena, 
similar to those of reactions (9), (12), and (13), will 
shift the apparent pKa in a manner of the above 

5 discussions. For example, the additives in drug 
formulations, such as surfactants, bile salts, 
phospholipids, ion-pair forming counterions, 
cyclodextrins or polymers may make the drug molecule 
appear more soluble. As long as such excipients do not 

10 undergo a change of charge state in the pH interval of 
interest (i.e., the excipients are effectively non- 
ionizable) , and the drug molecule is ionizable in the pH 
interval of interest, the difference between the 
apparent pKa and the true pKa will reveal the true 

15 aqueous solubility, as if the . excipient were not 
present - 

Table 1 summarizes some of the relationships 
developed between solubility, pKa, and pKa'^^. 



/ 
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TABLE 1 True Aqueous Solubility Detenained from pKa 
Shifts of Monoprotic Compounds 



lonizable 

Compound 
Type 


A = 


L.X Lie; 

aqueous 


Examples ( from 
Fig. 8 and 

XI ueira. uuic } 


acid 


A > 0 


log So"^^ - A 


diclofenac, 
furoseiuide, 

indomethacin, 
probenecid, 

naphthoic acid 


acid 


A < or « 0 


log So*^^ 


prostaglandin F2a ^ 


base 


A > or .= 0 


log So'^^^ 


phenazopyridine 


base 


A < 0 


log So''^^ + A 


amitriptyline, 
chlorpromazine, 
miconazole, 
terfenadine 



1680-1685, 



(AUC) DETEBMINED BY WEIGHTED 



AREA- UNDER-THE-CURVES 

REGRESSION METHOD ' 
10 The method utilized by the analytical device for 

determining concentrations of species by UV spectroscopy 

is based on a weighted regression analysis of whole 
-spectra. All the 65-130 measured absorbances from each 

microtitre plate well may be used in the analysis. The 
15 concentration of the reference species is known, and the 

object of the analysis is to assess the unknown 

concentration of the sample in the filtered solution, by 

appying Beer's law. 
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10 



15 



The calculated absorbance of the sample species at 
a particular wavelength i, ai can be related to the 

absorbance values of the reference species, ai , the 
absorbance of DMSO (scaled to the units of the DMSO 
concentration in the reference well), ai^°^^^\ and the 
absorbance of the sample solution ^^blank" (buffers and 
the plastic microtitre plate), ai^"''^*^''''' : 



ai 



^-^^ = po ai"''^ + pi ai"'^^'^^' + P2 ai^"'^''*"^' (14) 



where po, Pif and p2 are unknown parameters. In turn, 
the absorbance of the reference species is derived from, 
the observed reference spectrum, corrected for DMSO and 
blank: 



The ^reference' spectrum is that of the compound of 
known concentration, and the ^sample' spectrum is that 
20 of the same compound, but of unknown concentration. At 

. - DMSO (BEE) 2 - 

each • desired wavelength the value ot ai is 
independently established: the DMSO absorbance curve is 
characterized from measurements of solutions containing 
all components except for the sample species. The 

25 'blank' absorbances are measured before ea'ch assay, by 
sampling the microtitre plate filled only with the 
uniyersal buffer solution. So, ai ^ and ai 
are then known. 

The object of the analysis is to determine the 

30 three unknown parameters pof Pi, and p2. (The parameter 
Po is used directly in solubility determination; usually 
the other two parameters ' are determined to have near 
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unit values.) This is done by weighted linear 
regression analysis [P. R. Bevington, "'Data Reduction 
and Error Analysis for the Physical Sciences/' McGraw- 
Hill: New York, 1969, pp. 164-186], as outlined below, 
5 Define a column-vector of the unknown parameters, 

P = (Po/ Pi/ P2) . Define a row-vector, bi, at each 
wavelength i, consisting of the three independent . 
parameters ai^^ ^msom.)^ ^^msm^:^)^ 

10 bi = { ai^^ a^DMso(RBF)^ 3.B1.MSMIPLE) ) (16) 

so that equation (14) may be expressed in vector 
notation as 

15 ai -^*= = bi p (17) 

For all the measured wavelengths i (typically, i=0 to N- 
1, where N=number of wavelengths measured), equation 
(17) may be expressed in an matrix notation. 



20 



a ^^^^ = B p (18) 



Without the subscript, a^^-^^ denotes the column-vector of 
N rows (the number of discrete wavelengths measured) , 
25 and B denotes the N-rows x 3-columns matrix composed of 
the N individual bi vectors. 

; The parameters p are determined by finding a 
solution which minimizes the sum of weighted residuals 

30 ^ Q - Z i«o..N-i (ai^^^ - ai^^^^)2 / 

\ = a°^^'^ W a°^^ - p^ b' W a''''^ - a''^^''^ W/fi p (19)* 

+ p"^ b"^ W B p 
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where ai is the measured absorbance of the sample at 
wavelength i (dependent variable), is the estimated 
variance in the measured dependent variables ai**% T is 

5 the transpose operator, and W is the diagonal matrix of 
weights, each i'* element being the inverse variance, 
l/ai2. The model equation, ai""^", 'is a function of the 
three p parameters, as well as the independent 
variables, Xl (i^^ wavelength) and bi. 

10 The elements of p are determined by setting the 

derivative of Q with respect to p' to zero and solving 
for p. 

p = (B'^ W B)-^ (B^ W a"''^) (20) 



15 



20 



Note that po may be substituted into equations (5) and 
(7) as, 

Po = AUCs / AUCr (21) 



WEIGHTING SCHEME BASED ON PEAK SHAPE ANCMALIES 

The method utilized by the analytical device for 
determining concentrations of species by UV spectroscopy 

25 is based on a weighted regression analysis of whole 
spectra, as described in the last section. The 
technique is unique in how it assigns weights to the 
data, to make the assessment of concentrations more 
reliable, especially when samples are not entirely pure. 

30 The normal weighting scheme used in equation ^ (19) 

is constriiicted from the variances 
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jy2 (a) = + ( tJx da/d%) * 



(22) 



in the software of the analytical device, ac = 0.0002 
(absorbance units), the fixed contribution to the 
5 variance in the measured absorbance (experimentally 
determined by replicate baseline measurements, usually 
in the domain 450-500 nm) , and ctx = 0.2 nm (estimated 
error in wavelength, specified by the UV 
spectrophotometer manufacturer) . The dependent. 

10 variables, ai^S and the independent variables, ^i, are 
all assumed to be subject to error. The weighting 
scheme properly recognizes that measurements of 
absorbance on the steep sides of peaks are not as 
reliable as those near the peak top or from the baseline 

15 portions. 

However, suppose that the sample spectrum has an 
-impurity" absorption (or scattering due to dust 
particles) at a particular wavelength, .which is not 
found in the reference spectrum. If the impurity 
20 portion in the sample spectrum were not treated in any 
special way, then the procedure of the previous section 
would overestimate po, due to the added contribution of 
the impurity absorbance. 

Consider a sample solution containing a mixture of 
25 components (i.e., impurities). The observed spectrum 
would be a manifold of several contributing spectra. It 
is reasonable to assume that the solubilities of all the 
components are different. If a saturated solution is 
made of the mixture, some components will precipitate 
30 and some will not. The consequence of this is that the 
UV spectrum of the manifold may change in shape as a 
result of the precipitation. The shape-based weighting 
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scheme can detect such anomalies, minimizing their 
effect on the p-parameter determination. 

The calculation scheme assesses the peak maximum, 
a„aK'^^ and the wavelength, at which the maximum 

occurs for the reference spectra. For any given 
wavelength i, the following reference and sample ratios 
are calculated: 

r,^ = ai^^ / a^x'^^ (23) 



where a^^^'^^ is the sample absorbance at If the 

peak in the sample had the identical shape as the peak 
15 in the reference solution, then the two ratios, (23) and 
(24) would be the same. However, if it happens .that the 
impurity absorbance is at its maximum at the particular 
wavelength i, then the reference ratio would be greater 
than the sample ratio. 
20 Define a constant that is the ratio of the two above 

ratios : 

Si = 1 ri=^^= / ri^^ I (25) 



25 



30 



Furthermore, define the empirically-optimized scale 
factor kii 

ki = 10*®-' ' ^' for Si > 1 (26a) 

or 

^ IQ*6.3 (i/s - 1) for Si < 1 26b) 



- 34 - 

I 



wo 01/55698 



PCTAJSOl/02377 



now With this model, the normal weighting scheme in 
equation (19) is modified: 



20 



Cshape^ (a) = kM + { da/dX) - 3 



(27) 



Higher errors are ascribed to the spectral regions 
exhibiting impurities. Hence data from the regions of 
the impurity are made less important in the derivation 
o£ the parameter vector p, equation (20) . 

10 The way that the scale factor, k, is defined, it 

doesn't matter whether the impurity is in the sample or 
the reference spectrum: any differences in shape between 
the two spectra automatically down-weight the region of 
the differences in the regression equation (20), placing 

15 more emphasis on the spectral regions with the common 
shape . 

The following examples are intended to illustrate, 
but not limit, the invention disclosed herein. 

Example 1. Apparatus of Figs. 2 and 3: Solubility 
Determined by the Aqueous Dilution Method 

Fig. 2 summarizes the apparatus detailed below. 
The general aqueous dilution method used to determine 
25 solubility is described in Fig. 3a, and a specific 
embodiment is highlighted in Fig. 3b and detailed below. 

Data Coliectioji 

1. The analytical device places 2.5 . mL of 
30 universal buffer solution (25) , initially at pH 3,^" into 
the side arm tube of, a 'tall test tube 22 (-J-tube") 



wo 01/55698 



PCTAJSOl/02377 



10 



located in the pH titrator assembly (3). An 
alkalimetric pH titration is performed. 

2. Analysis of the titration data produces the pH 
values corresponding the titrant volumes 8, 16, 23, 30, 
37, 44, 52, 59, 65, 73, 81 nL; these correspond closely 
to pH settings 3.0 - 8.5, in increments of 0.5 pH units. 
These volumes will be used by the device below. 

3. The 96-well stock plate is placed in position" 
13 (S in Fig. 3b), furthest from operator on left rack 
on the robot table. Only one column of 8 sample wells 
will be used in this particular assay; for example, 
wells Al, Bl, HI (or A2, B2, H2, etc.) will 
contain the required compounds, 100 mM in DMSO in this 
example. This is the sample stock plate. 

15 4. A new microtitre plate (must be inert to DMSO - 

- NUNC 0.5mL polypropylene is satisfactory) is placed in 
position 14 (Sd in Fig. 3b) rear position of the rack on 
the right. This is the empty reference stock plate. 

5. A new 96-well deep plate (2.2 mL wells) is 
placed at position 15 (Dc in Fig. 3b), middle of the 
rack on the left on robot table. 

6. A new 96-well deep plate (2.2 mL wells) is 
placed at position 16 (Dd in Fig. 3b), middle of the rack 
on the right on robot table. 

7. A new 96-well Costar UV plate (or Greiner, but 
Costar is slightly better) is placed in position 17 (U) , 
nearest to operator on left the rack on robot table. 

8. A new rack of 20-fiL disposable (conductive) 
pipet tips is placed into position 10 and a new rack of 
200-jiL disposable (conductive) pipet tips is placed in 



20 



25 



30 



position 11 on the robot table. 
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9. A plastic waste bag is attached to the bottom 
of the waste-tip slide {12) . 

10. The system, solution bottle needs to have at 
least 500 mL of the universal buffer solution. 

5 11. A clean test tube in position 20, third test 

tube position nearest the operator of the test tube 
rack, is filled with 21 mL of DMSO. A clean test tube 
in position 21r second test tube position nearest the 
operator, is filled with 10 mL freshly prepared, low- 
10 carbonate, 0.5 M KOH. 

12. Four clean test tubes, each containing 10 mL 
distilled water, are placed into positions 18. 

•13. The robot loads 1000 universal buffer 

solution into each of the deep 96-well plate in position 
15 13 (DC in Fig. 3b) on the table. Then the robot draws 
0.5 M KOH from test tube in position 21 of the test tube 
rack. The robot proceeds to deposit 8 ^iL into the 8 
deep wells A2-H2, 16 jiL into deep wells A3-H3, 23 pi 
into deep wells A4-H4, 30 nL into deep wells A5-H5, 37 
20 nL into deep wells A6-H6, 44 ^L into deep wells A7-H7, 
52 ]iL into deep wells A8-H8, 59 jiL into deep wells A9- 
H9, 65 \iL into deep wells AIO-HIO, 73 pL into deep wells 
All-Hll, and 81 \il> into deep wells A12-H12. 

14. Robot loads an additional 1000 jiL universal- 
25 buffer solution (25) into each of the deep wells, 

followed by extensive mixing. 

15. Before any sample is added to 15 (Dc) , the 
robot transfers 150 jiL buffer solutions from the deep 
well to the UV, plate in position 17 (U) on the table (Dc 

30 — > U) . The operator is prompted to take UV spectra of 
the buffer-filled UV plate, using the spectrophotometer 
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(2) . Afterwards, the UV plate is rinsed with methanol- 
water using the plate washer (5) and returned to_ 
position 17 (U), covered with_a lid to keep dust out. 
This corresponds to the "blank" reading in equations 

(14) and (15) . 

16. TO each of the deep wells (containing 1850- 

1931 jiL pH buffer) , 5 jiL sample (27) is added, using the 
20 |il tips 13->15 (S -> DO) . The first sample pickup 
is discarded back into the stock plate, to ensure that 
the tip does not have an air gap at the opening and to 
ensure the dispenser stepper motor gears are free of 
slack, and that the inside surface of the tip is wetted 
before subsequent sample pickup. Vigorous 
-regurgitative" mixing follows. A protective cover is 
15 placed over 13 (S) plate after the transfers are 
complete . 

17. The top of the deep well plate 15 (Dc) is 
sealed (aluminum plate seals), and the plate is placed 
on the orbital shaker {6) for 4 hours. 

18. The robot transfers 190 ^lL DMSO to the empty 
reference stock plate 20->14. The robot transfers 10 
jiL samples: 13->14 (S -> Sd) . Solutions in Sd are 
thoroughly mixed. The sample stock plate 13 (S) may now 
be removed, protectively sealed, and stored in the 

25 refrigerator . 

19. The robot loads 1000 |iL universal buffer 
solution (25) into each well of the deep 96-well plate 
in position 16 (Dd) on the table. Then the robot draws. 
0.5 M KOH from test tube in position 21 of the test tube 
rack and proceeds to deposit 8 pL into the 8 deep wells 
A2-H2, 16 ^lL into deep wells A3-H3, 23 into deep 
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wells A4-H4, 30 nL into deep wells A5-H5, 37 nL into 
deep wells A6-H6, 44 into deep wells A7-H7, 52 \xL 
into deep wells A8-H8, 59 \iL into deep wells A9-H9, 65 
iiL into deep wells AIO-HIO, 73 |iL into deep wells All- 
5 Hll, and 81 nL into deep wells A12-H12. 

20. The robot loads an additional 1000 [iL 
universal buffer solution (25) into each of the deep 
wells, followed by extensive mixing. 

21. The robot discards 150 pL of each of the deep 
10 well solutions to waste. A second blank is not needed, 

but volumes in the deep well plate 16 (Dd) need to be 
the same those of IS (De) after step 13. 

22. To each of the deep wells, 15 \iL DMSO are 
added 20—>16. 

15 23. To each of the deep wells (containing 1850- 

1931 nL pH buffer), 5 [iL 20-fold diluted sample are 

added 14-->16 (Sd — > Dd) . Vigorous mixing follows, a 

critical step. 

24. As soon as mixing is complete in step 23, the 

20 robot transfers 150 nL of the reference solution to the 
UV plate 16->17 (Dd ->U). Filtration is not 
Necessary. The operator is prqn«>ted to take UV spectra 
of the UV plate, using 2. The deep well plate in 
•position 16 (Dd) is no longer used. The UV plate is 

25 thoroughly washed on the plate washer 5 (using 75% 
methanol - 25% water) and placed in position 17 (U) , 
with a protective dust cover. 

25. After the 4 hour shaking period elapses, the 
operator transfers the deep well plate from the orbital 

30 shaker to position 15 (Do) . 
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26. The operator places a used 96-well deep plate 
(2.2 itiL wells) in the bottom of the vacuum manifold, 
position 4 on the robot table. On top of that the 
operator places a new 96-well deep filter plate (2.2 mL 
wells) . The operator assembles the vacuum manifold 
accordingly. The robot transfers 20 precipitated 
solutions for the filtration step: 15->4. The 
operators is prompted to turn on the vacuum, to filter 
the solutions. This step pre-conditions the filters 
with sample, and is part of the so-called "double- 
filtration" step. 

27. The operator replaces the bottom (used) deep- 
well plate with a clean empty 96-well deep plate (2.2 mL 
wells) in the bottom of the vacuum manifold 4. On top 
of that the operator returns the sample-preconditioned 
96-well deep filter plate (2.2 mL wells). The operator 
assembles the vacuum manifold again. 

28. The robot transfers 250 [iL saturated sample 
solutions for the filtration step: 15->4. Afterwards, 
the deep well filter plate in position 15 (Dc) may be 
sealed and stored for future examination or re-use. 

29. The operators is prompted' to turn on the 
vacuum, to filter, the solutions, some of which may have- 
solids in them. The operator discards the top deep-well 
filter plate and exposes the bottom deep-well plate for 

robotic access. 

, 30. The operator removes cover on the UV plate at 
position 17 (U). The robot transfers 150 fxL to the UV 
plate 4— >17. The operator is prompted to take UV 
spectra of the UV plate. The deep well plate in 
position 4 may be discarded afterwards. If the plate 15 
{DC) is to be re-sampled. later, then the UV plate should 
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be thoroughly washed on the plate washer 5 (using 75%' 
methanol - 25% water) and placed in position 17 (U) , 
with a protective cover, for future use. Blank and 
reference spectra need not be repeated in such re- 
5 samplings . 
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Data Processing 

Fig. 9 shows the measured absorption spectra due to 
piroxicam. The shown absorption spectra had been 
corrected for the absorption of "blank" plate, and for 
5 the absorption due to DMSO (equation 14) . 

For ionizable compounds, spectra can change 
dramatically as a function of pH. Also, at a given pH, 
if a portion of the sample precipitates, the optical 
density (OD) will decrease, following Beer's law. 
10 Fig. 9 shows how spectra change with pH in the 

sample and reference solutions of piroxicam, a molecule 
which is effectively an ionizable acid if only the pH > 
3 is considered, with a pKa of 5.07 (Table 2). Both the 
sample spectra and the reference spectra are of 
15 excellent quality; however, one notes that the optical 
densities of the reference absorbance curves are 
relatively low (about OD 0.06 at the maximum), compared 
to those of the sample. If we had started out with 10- 
itiM DMSO stock solutions, instead of 100 mM, the CDs 
20 would not have exceeded 0.006, a useable value, but one 
nearing the limit of detection. To produce a reference 
solution free of precipitate, it is necessary to 
compromise on sensitivity. Example 1 conditions have 
been demonstrated to work well. However, if 10 mM DMSO 
25 stock solutions are used, different DMSO stock plate 
dilutions need to be considered (S ~> Sd) , along with 
different sampling volumes (S ~> Dc, Sd — > Dd) . 

Example 2 below will show how the sensitivity issue 
can be avoided by using cosol vents. 
30 As precipitation takes place to varying degrees in 

different pH values, the spectra of the sample solutions, 
change in optical densities. This can be clearly seen 
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in Fig. 9, where the sample spectra have the lowest OD 
values at pH 3.5 and systematically show higher OD 
values as pH is raised. The changing OD values indicate 
that solubility changes with pH. 
5 The reference solutions only changes as a result of 

the ionization of the compound, and not due to changes 
in concentration of the compound. Usually ionization 
does not produce large changes in the OD values as a 
function of pH. 

10 In the above example, the calculated reference 

solution concentrations in plate 16 (Dd in Fig. 3b or 35 
in Fig. 3a) are 

Cr = 10/ (190+10) X 5 / (2000+Vkoh -150+5+15) x 100 raM 
15 = 12.8 - 13.4 JIM ^28) 

where Vkoh, the volume of 0.5 M KOH added to adjust the 

pH, varies from 0 to 81 

The reference concentrations are indicated by the 
20 lower horizontal lines in Fig. 5. 

The liitLiting sample concentrations (assuming no 
precipitation had taken place in the deep-well plate 15 
(31 in Fig. 3a or Do in Fig. 3b) are calculated 
according to the formula 



25 



30 



Cs = 5 / (2000 + VKcm " 150 + 5) x 100 mM 

.• = 258 - 270 pM (assuming no precipitate forms) (29) 

The limiting concentration is represented by the 
horizontal lines in Fig. 7 and the upper horizontal 
line's in Fig. 5. 
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The next step requires that the spectra taken from 
190 to 500 nm (in increments of 2 nm) be analyzed, to 
determine the area-under-the-curve (AUC) values, using 
equation (20) incorporating the shape-based weighting 
5 scheme, equation (27), in the analysis. With plastic UV 
plates, usually only the data in the interval 240 - 500 
nm are actually used in the analysis. The optimized 
regression parameter is related R, the ratio of AUC 

values, equation (5) . 

If R = 1 and the .reference solution is not 
supersaturated, then .one cannot report the solubility of 
the compound at the particular pH, since both the sample, 
and the reference solutions have precipitation in them. 
This is indicated hypothetically in Fig. 5A for pH < 
PHk, and in Fig. 6A by the low-pH region flat portion of 
the curve. 

This situation is seldom seen in practice, since 
the quick handling of the reference plate and the higher 
amount of DMSO used in the reference plate compared to 
the sample plate almost always ensures that the 
sparingly soluble compounds are supersaturated. This is 
made more probable by the use of 1.06% v/v DMSO in the 
reference plate, compared to the usual 0.26% v/v DMSO in 
the sample deep well plate (in the . latest variant of the 
protocol). On the other hand, the sample plate is 
allowed to sit for at least 4 hours, ample time in most 
cases for precipitation to start. The longer the wait 
time, the more "crystalline" the precipitated solid 
becomes . 

30 An earlier version of the protocol used 0.5% v/v 

DMSO in both the reference and the sample deep well 
plates. The latest protocol (1.06% v/v DMSO in 
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reference and 0.26% in sample) enhances the chances of 
supersaturation conditions in the reference deep well 
plates . 

If R = 10 / (190+10) = 0.05, then one cannot report 
5 the solubility of the compound at the particular pH, 
since both the sample and the reference solutions have 
no precipitate. This condition is common, and is 
indicated for all pH in Fig. 7A, for pH < 8 in Fig. 7B, 
for pH < 7 in Fig. 7C, for pH > 6 in Fig. 7D, and for pH 
10 > 6 in Fig. 7D. This corresponds to the hypothetical 
conditions for pH > pHs in Fig. 5A and Fig. 5B, and flat 
regions at the lowest parts of the curves in Fig. 6. To 
remedy this, one needs to use more concentrated DMSO 
stock solutions. This is not a problem in practice, 
15 since many discovery applications are mainly concerned 
with Cs < 50 ^g/mL. We demonstrated the aqueous, 

dilution method with Cs of 100 and 200 fig/mL (Fig. 7). 

The R vs. pH relationships (Fig. 6) can be 
automatically assessed by the analytical device, and 
20 when the R is less than 1 but greater than 0.05, 
solubility is calculated according to equation (6) . The 
useful pH domain for solubility determination are 
indicated in Fig. 5 as pHr < pH < pHs- In Fig. 7B, 
solubilities of amitriptyline can be determined for pH > 
25 8; in Fig. 7C, solubilities can be determined for pH > 
7 , etc . 

The minimum possible solubility refers to the 
intrinsic solubility, and such So values often need to 
be extrapolated from the S vs. pH curve, taking its. 
30 theoretical expected shape into account. The apparent 
intrinsic solubilities, So"^ determined in this way are 
listed in Table 2, column 3, for the compounds used in 
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the example. All the So^^' values reported in Table 2 
were determined in the presence of 0.5% v/v DMSO, except 
for phenazopyridine, where 0.26% was used. 

5 EesuJts of Aqueous Dilution J^lication 

The most insoluble compounds were found to be 
terfenadine, miconazole, indomethacin, chlorpromazine, 
and piroxicam. We saw no evidence of precipitation with 
amiloride and nortriptyline, which suggests their 

10 solubilities are greater than 208 jig/mL and 98 fig/mL, 
respectively, in 0.5% v/v DMSO. Although we did see a 
few crystals in the propranolol well at pH 9.5, we were 
not able to derive the solubility with any confidence. 
We estimate that the solubility of propranolol is about 

15 100 ng/mL in the presence of 0.5% v/v DMSO. 

Example 2. Apparatus of Figs. 2, 4a, and 4b: Solubility 
Determined by the Cosolvent Method 

Fig. 2 summarizes the apparatus detailed below. 
20 The general cosolvent method used to determine 
solubility is described in Fig. 4a, and a specific 
embodiment is highlighted in Fig. 4b and detailed below. 
In the variation utilizing cosolvent, all spectra were 
collected of solutions containing 50% v/v 1-propanol, 
but in other respects the methods are very similar. 
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30 



Data Collection 

1. The analytical device places 2.5 mL of 
universal buffer solution (25), initially at pH 3, into 
the side arm tube of a tall test tube 22 («J-tub6") 
located in the pH titrator assembly (3) . - An 
alkalimetric pH titration is performed. 
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2. Analysis of the titration data produces the pH 
values corresponding the titrant volumes 8, 16, 22, 30, 
37, 44, 52, 59, 65, 73, 81 nL; these will correspond 
closely to pH settings 3.0 - 8.5, in increments of 0.5 

5 pH units. These volumes will be used by the device 
below. 

3. The 96-well stock plate is placed in position 
13 (S in Fig. 4b), furthest from operator on left rack 
on the robot table. Only one column of 8 sample wells 

10 will be used in this particular assay; for example, 
wells A1, B1, HI (or A2, B2, H2, etc.) will 

contain the required compounds, 100 mM in DMSO' in this 
example. This is the sample stock plate. 

4. A new microtitre plate (must be inert to DMSO - 
15 - NUNC 0.5mL polypropylene is satisfactory) is placed in 

position 14 (Sd in Fig. 4b) rear position of the rack on 
the right. This is the empty reference stock plate. 

5. A new 96-well deep plate (2.2 mL wells) is 
placed at position 15 (Dc in Fig. 4b), middle of the 

20 rack on the left on robot table. 

6. Test tube in position 19 [42 in Fig. 4a) is 
filled with 20-21 mL spectroscopic-grade 1-propanol. 

7. A new 96-well Greiner UV plate (Greiner is 
better than Costar in- the cosolvent method) is placed in 

25 position 17 (U) , nearest to operator on left the rack on 
robot table. 

.- 8. A new rack of 20 disposable (conductive) 
pipet tips is placed into position 10 and a new rack of 
200 ]iL disposable (conductive) pipet tips is placed in 
30 position 11 on the robot table. 

9. A plastic waste bag is attached to the bottom 
of the waste-tip slide (12) . 
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10. The system solution bottle needs to have at. 
least 500 mL of the universal buffer solution. 

11. A clean test tube in position- 20, third test 
tube position nearest the operator of the test tube 

5 rack, is filled with 21 mL of DMSO. A clean test tube 
in position 21, second test tube position nearest the 
operator, is filled with 10 mL freshly prepared, low- 
carbonate, 0.5 M KOH. 

12. Four clean test tubes, each containing 10 mL 
10 distilled water, are placed into positions 18. 

13. The robot loads 1000 jiL universal buffer 
solution into each of the deep 96-well plate in position 
13 (DC) on the table. Then the robot draws 0.5 M KOH 
from test tube in position 21 of the test tube rack. The. 

15 robot proceeds to deposit 8 ^L into the 8 deep wells 
A2-H2, 16 jiL into deep wells A3-H3, 23 tiL into deep 
wells A4-H4, 30 JiL into deep wells A5-H5, 37 into 
deep wells A6-H6, 44 [iL into deep wells A7-H7, 52 
into deep wells A8-H8, 59 pi. into deep wells A9-H9, 65 
20 ^tL into deep wells AlO-HlO, 73 pL into deep wells All- 
Hll, and 81 |iL into deep wells A12-H12. 

14. Robot loads an additional 1000 ^lL universal 
buffer solution (25) into each of the deep wells, 
followed by extensive mixing. 
25 15. Before any sample is added to 15 (Do), the 

robot transfers 73 ^ buffer solutions from the deep 
well to the OV plate in position 17 (U) on the table (Dc 
-> u) ; the robot transfers 72 pL of 1-propanol from 15. 
The robot mixes the solvents in each well. 
30 16.' The operator is prompted to take UV spectra of 

the cosolvent-buffer-filled UV plate, using the 
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spectrophotometer (2). The solution is SAVED and the 
plate is returned to position 17 (U) , with its top 
tightly covered, to keep the dust out and evaporation .. 
minimized. 

5 17. The robot transfers 190 pL DMSO to the empty- 

reference stock plate 20->14. (1-propanol is not 
suitable at this point because of its volatility.) The 
robot transfers 10 HL samples using the 20 pi. tips: 
13— >14 (S — > Sd) . Solutions in 14 (Sd) are thoroughly 

10 mixed. 

18. As soon as the solutions are mixed in step 17, 
the robot transfers 5 ^lL from the diluted sample stock 
solution (using the 20 |iL tips) to the blank solutions 
whose spectra had just been read in step 16: 14— >17 (Sd 
15 -> U) . There are now 150 \iL in each well of the Greiner 
plastic UV plate. The diluted sample stock plate 14 
(Sd) may now be removed, protectively sealed, and stored 
in the refrigerator. 

19. The operator reads the UV of the plate with 
20 the spectrophotometer 2. This becomes the reference 

plate reading. Discard the solutions and wash plastic UV 
plate thoroughly on the plate washer 5 (75% methanol - 
25% water). Place empty clean UV plastic plate in 
position 17 (U). Temporarily cover with lid to keep 

25 dust out. 

20. To each of the deep wells at 15 (Dc) 

(containing 1927-2008 nL pH-adjusted universal buffer), 
5 pL sample is added, using the 20 tips 13-->1S 

(S — > DC) . The first sample pickup is discarded back 
30 into the stock plate, to ' ensure that the tip does not 
have an air gap at the opening and to ensure the 
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dispenser stepper motor gears are free of slack, and 
that the inner walls of the tip are wetted with sample. 
Vigorous., -regurgitative" mixing follows. The sample 
stock plate 13 (S) may now be removed, protectively 
sealed, and stored in the refrigerator. 

21. seal the top of the deep well plate (aluminum 
plate seals), and place on the orbital shaker {6) for 4 
hours . 

22. The operator places a used 96-well deep plate" 
(2.2 raL wells) in the bottom of the vacuum manifold (4). 
On top of that the operator places a clean 96-well deep 
filter plate (2.2 mL wells) . The operator assembles the 
vacuum manifold accordingly. The robot transfers 20 
precipitated solutions for the filtration step: 15— >4. 
The operators is prompted to turn on the vacuum, to 
filter the solutions. This step, a part of the "double, 
filtration" operation, pre-conditions the filters with 
sample . 

23. The operator replaces the bottom (used) deep- 
well plate with a new 96-well deep plate (2.2 mL wells) 
in the bottom of the vacuum manifold. On top of that 
the operator returns the sample-preconditioned 96-well- 
deep filter plate (2.2 mL wells). Operator assembles the 
vacuum manifold again. 

24. The robot transfers 250 jxL precipitated sample 
solutions for the filtration step: 15->4. Afterwards, 
the- deep well filter plate in position 15 (Dc) may be 
stored for future examination or re-use. 

25. The operators is prompted to turn on the 
30 vacuum, to filter the solutions, some of which may have 

solids in them. The operator discards the top deep-well 
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filter plate and exposes the bottom deep-well plate for 
robotic access. 

26. The operator removes cover on the UV plate at 
position 17 (U) . The robot transfers 75 \iL of the 

5 filtered sample solution to the UV plate: 4— >17. 
Immediately following that, the robot transfers 75 jiL 1- 
propanol to the UV plate: 15~>17. The solutions are 
thoroughly mixed. 

27. The operator is prompted to take UV spectra of 

» 

10 the UV plate (17~>2) . The UV plate and the deep well 
plate in position 4 are discarded afterwards. 
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Data Processing 

Fig. 10a shows the measured absorption spectra of 
miconazole, and Fig. 10b shows those of phenazopyridine . 
The shown absorption spectra had been corrected for the 
5 absorption of the "blank" plate, and for the absorption 
due to DMSO (equation 14) . 

For ionizable compounds, spectra can change 
dramatically in shape as a function of pH, as in the 
reference spectra in Fig. 10b, but this isn't always so, 
10 as illustrated by the spectra of miconazole in Fig. 10a. 

The reference solutions only changes as a result of 
the ionization of the compound, and not due to changes 
in concentration of the compound. 

Also, at a given pH, if a portion of the sample 
15 precipitates, the optical density (OD) will decrease, 
following Beer's law, which is what all the sample 
spectra in Figs. 9 and 10 show. 

Fig. 10a shows how spectra change with pH in the 
sample solutions of miconazole, an ionizable base with a 
20 pKa 6.07 (Table 2). Both the sample spectra and the 
reference spectra are of excellent quality, with the 
optical densities of the reference absorbance curves 
being comparable to those of the sample. This is one 
of the advantages of the cosolvent method over the 
25 aqueous dilution method, in that sensitivity is similar 
in the sample and reference solutions. One of the 
disadvantages of the cosolvent method is that of the 
increased evaporation of the cosolvent-containing 
solution. Nevertheless, example 2 conditions have been 
30 demonstrated to work well. However, if 10 itiM DMSO stock' 
solutions are used, different DMSO stock plate dilutions 
need to be considered (S — > Sd) . 
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■ As precipitation takes place to varying degrees in 
different pH values, the spectra of the sample solutions 
change in optical densities. This can be clearly seen 
in Fig, 10a for the sample spectra, where the sample 

5 spectra have the lowest OD values at pH 9.0 and 
systematically show higher OD values as pH is lowered, a 
pattern consistent with that of an ionizable base. The 
changing OD values indicate that solubility changes with 
pH Fig. 10b is that of an ionizable acid, and so the 

10 pH - OD trend in the sample spectra follow the reverse 
order. 

in the above example, the calculated reference 
solution concentrations, after the cosolvent had been 
added in step 18 {48 in Fig. 4a) above, are 

15 

^co3 ^ (190+10) X 5 / (145+5) x 100 mM 
=167 MM 

Note that this is nearly ten times higher than the value 
20 in the aqueous dilution method (equation 28) . 

The limiting sample concentrations (assuming no 
precipitation had taken place .in the deep-well plate 15 
(DC) and in solution 43 in Fig. 4a) are calculated 
according to the formula 

25 

^cos ^ 75/150 X 5 / (2000 + Vkoh - 73 + 5) x 100 mM 

= 248 - 259 MM (assuming no precipitate forms) (31) 

The next step requires that the spectra taken from 190 
30 to 500 nm (in increments of 2 nm) be analyzed, to 
determine/ the area-under-the-curve (AUG) values, us.ng 
equation (20) incorporating the shape-based weighting 
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scheme, equation (27), in the analysis. With plastic DV 
plates, usually only the data in the interval 240 - 500 
nm are actually used in the analysis. The needed 
regression parameter is R'^-, the ratio of AUG values, 

5 equation (7). 

The same considerations need to be made with 
regards to the limiting cases of the R=°= ratios, as was 
discussed in the aqueous dilution method. 

The R = 1 situation (precipitate in both 43 and 48 
10 in Fig. 4a) is very rarely seen in practice, since the 
quick handling of the .reference plate and the higher 
amount of DMSO used in the reference plate compared to 
the sample plate almost always ensures that the 
sparingly soluble compounds are supersaturated if not 
15 dissolved. This is made more probable by the use of 
3.33% v/v DMSO in the reference plate, compared to the 
0.26% v/v DMSO in the sample deep well plate, and by the 
added cosolvent, 1-propanol at 50% v/v. In effect, 
there are fewer restrictions in the cosolvent method. 
20 volatility of the cosolvent is the main disadvantage, as 
said before. 

The solubility is calculated according to equation 
(8). The results for phenazopyridine stated in Table 2. 
were determined by the cosolvent method. 

Example 3. Aqueous Solubility Determined From the Shifts 
in ...pKa in the Presence of Non-ionizable Chemicals that 
Distort the Solubility-pH Profiles. 

The apparent solubilities determined by equations 
30 (6) and (8) are plotted in Fig. 7 as the points below 
the horizontal lines. Those points are replotted in 
Fig. 8 as log S*^^ vs. pH. The solid curves in Fig. 8 
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are the results of fitting the apparent solubilities to 
standard equations described in the Background section/ 
From that analysis, it was possible to obtain pKa'^^^ and 
So"" values. The latter apparent intrinsic solubility 
5 constants are listed in Table 2. 

As a visual check, the apparent ionization 
constants may be estimated as the pH values at the 
intersections of the horizontal and diagonal asymptotes 

of the solid curves. 
10- Since the pKaS of the studied compounds are known 

(Table 2), it was .possible to calculate the pKa shifts 
(A values in Table 1). These shifts were used to 
calculate the corrected aqueous intrinsic solubilities. 
So, also listed in Table 2. 
15 As one can see in all plots in Fig. 8," 

DMSO/aggregation has a significant influence on the 
intrinsic solubility even in 0.5% v/v DMSO. Most of the 
pKaS were shifted by up to one log unit. As far as we 
could find, this had not been recognized before in DMSO- 
20 containing solutions. Nor had a method been put forward 
that could correct for this effect (cf.. Table 1). Table 
2 summarizes the results of the corrective calculations, 
and compares the results to literature values. The 
■ comparison is especially good for terfenadine, 
25 probenecid, miconazole, and furosemide. 

•a. Chlorpromazine 

The ionizable base chlorpromazine has a pKa of 9.24. 
(Table 2) . From the solid curve in Fig. 8, the apparent 
30 pKa is 8.48, which is a pK^ shift of -0.76,, consistent in 
sign with the expected behavior of bases.. According to 
equation (11) , 
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10 



or 



15 



log So = log 19.4 ng/mL - 0-76 = 0.528 
So = 10*°-"" = 3.4 jig/mL 



(32) 



(33) 



The horizontal asymptote in the dashed curve 
corresponding to chlorpromazine in Fig. 8 is 0.528 (log 
scale), corresponding to 3.4 ng/mL. This value is in 
excellent agreement with the value determined with the 
pSOL instrument, and. is in reasonable agreement with the 
shake-flask value (Table 2). We cannot say with 
certainty that the pK^ shift is due solely to the effect 
of DMSO, or due solely to aggregation, since 
chlorpromazine is known to be surface-active. However, 
the method of our invention does not require the 
knowledge of the source of the effect. The method 
simply corrects solubilities according to the scheme xn 
Table 1. 
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b. Piroxicam 

The ionizable acid piroxicai; (considering pH > 3) 
has a pKa of 5.07 (Table 2). From_ the solid curve in 
Fig. 8, the apparent pKa is 6.05, which is a pK. shift of 
5 +0.98, consistent in sign with the expected behavior of 
acids. According to equation (10), 

log So - log 10.5 (ig/ioL - 0.98 = 0.041 (34) 

or 

10 So = 10"^-"' = 1.1 1^9/°^ ^^^^ 

The horizontal asymptote in the dashed curve 
corresponding to piroxicam in Fig. 8 is 0.041 (log 
scale), corresponding to 1.1 ^g/mL. This value is in 
15 poor agreement with the values determined by the shake- 
flask method (Table 2), 9.1 and 8-16 ^ig/mL. Both of the 
latter values were determined in the absence of DMSO. 
The inspection of the log S vs. pH data corresponding to 
the former value indicated a positive pKa shift. 
20 Application of the correction yielded the value 3.3 
Hg/mL, a considerable improvement. The 8-16 jig/mL. 
reported value did not reveal detailed log S vs. pH 
data. However, the authors did indicate that the pKa 
was deduced from the data as 5.63. This would suggest 
25 the expected positive shift, and the corrected intrinsic 
solubility is thus 2.2- - 4.4 fxg/mL, in very good 
agreement with the other 1.1 and 3.3 |ig/mL values. 
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TABLE 2 Intrinsic Solubility, So, Corrected for the 



compound 




(^xg/mL) 


corrected 
So 
([xg/mL) 


pSOL 

So 
{jig/mL) 


Shake-Flask 

So 

(lig/mL) 


aiaitriptyline 


9.45 ^ 


56.9 




2.0 ^ 


2.0 ' 


chlorpromazme 


9.24 ^ 


19.4 


3.4 


5 ^ 


0.1 ^ 


diclofenac 


3.99 " 1 


22.6 


3.8 


w • o 


0.6 " 


furosendde 


10.63, 
3.52 ^ 


29.8 


2.9 


5.9 " 


12.0 ° 
(2.9 ") 


griseofulvin 


non- 
ioniz- 1 

able 


37.6 


20.2 




9 » 


indomethacin 


4.42 " 


,7.2 


4.1 


2.0 ^ 


"2.0 % 1 ' 


miconazole 


6.07 ^ 


11.1 


1.6 


0.7 




2 -naphthoic 
acid 


"4.16 ' 


33.3 


20.2 




22.4 ' 


phenazo- 
pyridine 


■ 5.15 * 


12.2 


"1272 


■ 14.3 




piroxicam 


■ 5.. 07, 
2,33 ^ 


10.5 


"T7I 




(3.3 =), 

8-16^ 

(2.2-4,4°) 


probenecid 


3.01 " 


4.6 


0.7 


-0.6* 




terfenadine 


" 9.53 " 




0.1 jO-1 

n lir-l-nysRrin . C.A.S 


- 

. Johansson,. 



10 



15 



^ M A Strafford, A. Avaeei, r. ru-^^^-ww — , _ . 

K Luttoan, C.R. Brownell, R. Lyon, M. Assoc. Pbana. Sci. 
Ann. Mtng. 2000, poster presentation. 

" A. Avdeef, CM. Berger, C. Brownell, Pham. Res. -suuw, 

o'clr^ected for aggregate formation: ^^^l^'^^^.^f^'^ soa. 
<« J. Huskonen, M. Salo, J. Taskinen, J. Chem. Int. Conp. i>oc. 

"I'h l^lJow^kyf'R.-M. Dannenfelser (Eds.), AQOASOL dATAbASE 
of Aqueous S^lS-ility, 5th Ed., 1998, College of Pharmacy, 
Univ. of Arizona, Tucson, AZ 85721. 

^ pION, unpublished data. cteila. J. 

9KG. Mooney, M.A. Mintun, K.J. Hutffliestein, V.J. Stella, 

Pharm. Sci, 1981, 70, 13-22. 
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^ A. Avdaef, K.J. Box, Sirius Technical Application Notes, 
vol. 2. "96, p. 110^ correspondence, 2000. 

■5 ' Sci! 2000, eS, 268-274 (24 hr) . 

c . Phenazopyridine 

Phenazopyridine is an interesting molecule. It is 
an ionizable base, but it has positive pK. shifts. The 
10 true pKa is 5.15, but in saturated solutions, containxng 
0 26% v/v DMSO, the apparent pK. is 5.83. We believe 
that this is a case where positively charged species 
either interact with -DMSO "or form aggregates. Bexng 
Charged species, the effect is to shift the diagonal 
15 segment of the curve in a horizontal direction xn the 
log S vs. PH plot, as indicated by the dashed Une in 
Fig. 8. This means that the intrinsic solubility 
constant is not affected (Table 1) . 

we found a similar example in the literature for 
20 prostaglandin F2a (Table 1) . The molecule is an acid 
but Shows a negative pK. shift. The authors proposed 
that micelles formed to cause the shift. No effort was 
n.ade to characterize the apparent pKa. As in the case 
of phenazopyridine, the intrinsic solubility is not 
25 affected by the aggregation of negatively-charged 
species - 



/ 



- 59 - 



wo 01/55698 



PCT/USOl/02377 



10 



Kx»ple 4. Shape-fca«d «alghti.g Scheme Spiled to 
Solubility Determination 

T= 4r, ir^n 11. Shown are 
. Let's consider the example in Fig. ix- 

the spectra of a sample and a reference solution of 
piroxicam. The peak maximum occurs at 359 nm. Consider 
the point in the spectra at 285 nm, where there is a 
p.onounced difference in shapes between the reference- 
spectru. and the sa^^le spectrum. The —1 error at 
285 nm is calculated by equation (23) to be 0.0007 
absorbance unit. For this example, 

^^BKF ^ a^HEF / a^BE^ = 0.06/0.104 = 0.577 (36) 
^^SAMPLB » a.sAMPi^ / a^'«^ - 0.50/1.025 = 0.488 (37) 

Si = 1 ri""'''^ / I = O-'*^^ ^ ^'^"^"^ " ^^^^ 

= 10^«-^ ( x/s - 1 ) = 13 

, so, the Shape-biased error is a assigned to be 13 times, 
higher in the region of shape discordance, namely, 
0.008. Hence, the region near 285 nm is down-weighted 
by a factor of 13^ or 169. In effect, when peak shape 
differences are encountered, the data are automatically 

5 taken out of consideration in the determination of the 
p-vector in the weighted linear regression analysis. 
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Example 5. Shape-based Weighting Scheme Implied to 
Permeability Determination 

The shape-based weighting scheme we developed is 
general and can be applied to any spectroscopic 
determination -of concentrations using multi-wavelength 
data we have applied the UV-based shape analysis 
weighting scheme to the determination of permeability 
constants. This scheme is also applicable to any 
concentration determination using spectral data. 

EKample 6. Improved Aqueous Universal Buffer Formulatioxr 
possessing Low-UV Components and Being Free of Phosphate 
and Citrate 

in this example and in Example 7 below, two new 
improved universal buffers are described for high- 
throughput solubility applications, and other 
applications where concentration measurements are made 
at controlled pH using spectrophotometric methods r' 
including permeability measurements. We sought to 
discover an aqueous and a cosolvent-water universal, 
buffer which would make solubility measurements more 
reliable. We wanted buffers that would not interact 
with the molecules being studied, which meant that the 
background salt concentrations needed to be low. For. 
this reason, we had to stay away from phosphate as a 
buffer component due to its strong tendency to cause 
salts of positively-charged drug -substances to 
precipitate [W.H. Streng, S.K. Hsi, P.E. Helms, H.G.H. 
Tan, J. Pharm. Sci. 1984, 73, 1679-1684]. We wanted to 
30 avoid ion-pair interactions between buffer components 
and the sample, so we wanted to stay away from 
^ lipophilic buffers. Since we wanted our detection 
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system to obtain as much signal due to aa»ple and as. 
little due to background buffers, we wanted the W 
absorption due to the buffer to be as low as possible 
Citric acid buffer component would not have been a good 
choice from this consideration. We wanted the pH vs. 
volume of alkaline titrant relationship to be as linear 
as possible, so that we could easily predict how much of 
a pH change a given sample may impart on the buffer 

system • . 
None of the conononly known universal buffers [D.D. 

Perrin, B. Dexnpsey, Buffers for pH and Metal Ion 
control. Chapman and Hall, London, 1974; A. Avdeef, J- J. 
Bucher, Anal. Chen.. 1978,50,2137-2142] fit the desired 
profile, we found that the zwitterionic Good's buffers 
15 [Good, Winget, Winter, Connolly, Izawa, Singh 
Biocteraistry 1966,5,467] best fit our desxred 
characteristics. Several buffers were tested as possible 
candidates, but one set was discovered to be 
outstanding . 

20 The best candidates for the aqueous universal 

buffer were discovered to be the con^ination: acetic 
■ acid, MES(4-morpholine-ethanesulfonic acid), HEPES(4-(2- 
hydroxyethyl) -l-piperazineethanesulf onic acid) , and 
boric acid. MES and boric acid had virtually no UV 
25 absorption above 210 nm. HEPES and acetic acid had pK^s 
in the right positions to make them attractive choices, 
even though their UV absorptivities were not the lowest. 
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Each component was initially chosen to be 5 mM, to 
form an equimolar mixture, but better ratios were 
discovered afterwards. 
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a TO achieve a starting pH 3.0, lb.26 mL ^..M hux 

were added per liter, ^°^^^^P°^^,\"5.f J^^^^^^ 
HCl. I = average ionic strength m the titration rrom 

pH 3 - 10. 

Table 3 shows the best component concentrations, 
along with other properties of the . components. 
Obviously, other forms of the buffering materials (e.g., 
sodium acetate, anhydrous) could be utilized. Figure 12 
shows the pH-volume curve, where 2.00 xtCL of aqueous 
universal buffer solution were titrated with 0.5 M KOH. 
The^ linear regression of pH vs. titrant volume produced 
the equation: pH = 3.10 + 72.08 Vkoh. with the 
extraordinary standard deviation of 0.02 pH units. 



- 63 - 



PCTAJSOl/02377 

WO 01/55698 

EKample 7. Improved Aqueous -Cos o Ivan t Universal Buffer 
Formulation Possessing Low-W Components and Being Free 
of Phosphate. .and Citrate 

Finding suitable candidates for a cosolvent 
5 universal buffer solution was problematic, since 
cosolvents change the pK^s of the buffer components, and 
it is difficult to predict accurately to what extent the 
pKaS would be different. The best candidates for the 
cosolvent buffer (Table 4), where the cosolvent 
10 acetonitrile was mixed with water in 1:1 proportion 
were discovered to be glycolic acid, MES, HEPES, and 
taurine (2-aminoethanesulfonic acid) . 

Figure 13 shows the pH-volume curve, where 2.00 mL 
of acetonitrile-water universal buffer solution were 
15 titrated with 0.5 M KOH. The linear regression of pH 
vs. titrant volume produced the equation: pH = 3.01 + 
75.59 Vkoh, with a good standard deviation of 0.06 pH 
units . 
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Table 4 Con^osition of tHe 1:1 Acetonitrile-Water 
Universal Buffer 




: r~::rT;:Tn:irTo S 76 mL of 0.5M HCl were 

:i:J'':::ut%Toor^'^^^^ to o'.0O3366 « adde^ HCI 
?^?SaUv 500 io. aoetonitrile wa= added to the weighed 
i™en« in a 1 L vol.metTlo flask, to which were - 
Sed 6 76 S ?.5 M HCl. Distilled water was then added 
10 to make 1.00 L solution. 

variations of the above "methods and analytical 
devices will be apparent to those with ordinary skill in 
the field. 
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CLAIMS 

We claim: 

5 1 A method for determining the solubility of a 
compound, said compound having limited solubility in 
water and in a buffer, comprising: 

a preparing a saturated solution of saxd 
• compound, after which any insoluble compound is filtered 

b preparing a reference solution of said 
compound, under conditions avoiding or suppressing 

precipitation^ 

c. measuring a spectrophotometric property of 
15 said saturated and reference solutions and said buffer, 

wherein said buffer serves as a blank, and 

d. analyzing said spectrophotometric property to 
determine the solubility of said compound. 

20 2 The method of claim 1, wherein said, 
spectrophotometric property is selected from the group 
consisting of UV range spectrophotometry, visible range 
spectrophotometry, colorimetry, light scattering 
detection, polarimetry, optical rotation, fluorimetry 
25 and circular dichroism detection. 

3. •• The method of claim 1 wherein: 

a. said reference solution is prepared by 
dissolving a known quantity of said coit«.ound or stock 
30 solution containing said compound in buffer to produce a 
solution /free of precipitate, and 
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b said saturated solution is prepared by 
dissolving a known quantity of said compound or stoclc 
solution containing said compound in buffer to produce a 
solution with precipitate dispersed therein. 

4 The method of claim 1 wherein said compound is 
presented as a stock solution in DMSO and said 
spectrophotometric property is UV range absorbance. 

10 5. The method of claim 4 wherein: 

a. said reference solution is prepared by 
diluting said stock solution with DMSO and buffer, and 

b said saturated solution is prepared by 
diluting said stock solution with buffer, after which 
15 any insoluble compound is filtered off. 

6. The method of claim 5, wherein: 

a. said reference solution is prepared by 
diluting one aliquot of said stock solution with 19, 
20 aliquots of DMSO; and one aliquot of said diluted, 
solution is further diluted with 3 aliquots of DMSO and 
400 aliquots of buffer, 

b said saturated solution is prepared by 
diluting 1 aliquot of said stock solution with 400 
25 aliquots of buffer, allowing said solution to 
equilibrate for several hours, after which any insoluble 
compound is filtered off, and 

c. said analysis requires that (1) the absorbance 
of said reference solution be compared to (2) the 
30 absorbance of said saturated solution at one or more 
wavelengths, wherein^ all of said solutions are adjusted 
to the same pH. 
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1 The method of claim 4, wherein 

a. said blank is prepared by diluting an aqueous 

buffer with cosolvent, 
5 b said reference solution is prepared by 

diluting stock solution with DMSO, buffer, and 
cosolvent, 

c said saturated solution is prepared by 
diluting stock solution with buffer, filtering off any 
10 insoluble compound, and further diluting said solution 
with cosolvent- 

8 The method of claim 7 wherein said cosolvent is 
selected from the group consisting of acetonitrile, 
15 methanol, ethanol, iso-propanol, 1-propanol, 1,4- 
dioxane, dimethylformamide, acetone, ethylene glycol, 
propylene glycol, polyethylene glycol 400, 
tetrahydrofuran, DMSO and mixtures thereof. 

20 9. The method of claim 8, wherein said cosolvent is 1- 
propanol . 

10. The method of claim 9, wherein: 

a. said blank is prepared by diluting 1 aliquot 
25 of said buffer with 1 aliquot of a cosolvent, 

b said reference solution is prepared by 
diluting 1 aliquot of stock solution with 19 aliquots of 
DMSO; and 1 aliquot of said diluted solution is further 
diluted by 15 aliquots of said buffer and 15 aliquots of 

30 cosolvent; 

c. said saturated solution., is prepared by 
diluting 1 aliquot of stock solution with 400 aliquots" 
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of buffer, allowing said solution to equilibrate for 
several hours, after which any insoluble co-pound xs 
filtered off; and 1 aliquot of said resulting solution 
is further diluted by I aliquot of cosolvent; 
5 d. said analysis requires that (1) the absorbance of 
said reference solution be compared to (2) the 
absorbance of said saturated solution at one or more 
wavelengths, wherein, all of said solutions are adjusted 
to the same pH. 



10 



11. The method of claim. 1, wherein 

a. said solubility is intrinsic solubility, 

b. 



said compound is ionizable and anomalies are 
present in solution, and wherein said analysis requires 
15 that the log of solubility of said saturated solution be 
plotted against pH to produce a solubility-pH curve, 

c the displacement of the pK. caused by said 
anomaly can be determined by independently determining 
the pKa of said compound, and 
20 d. said displacement of pK, can be used to 

^v,<. effect of the anomaly and determine 
removing the eirecx; oi. v-i - 

intrinsic solubility of said compound. 

12 The method of claim 11, wherein said anomaly is 
25 caused by one or more factors selected from the group 

consisting of: 

a. self-association effects, 

" b" presence of bile acids or other surfactants, 

c. presence of cyclodextrins, 

30 d. presence of ion-pair forming counterions, 

e. presence of non-ionizable polymers, 

f . presence of phospholipids, and 
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g. presence of DMSO. 

13. The method of claim 12,. -wherein said anomaly is 
caused by the presence of DMSO, wherein: 
5 a said pKa of said confound in buffer free of 

anomaly can be determined by using known analytical 
techniques , 

b. the pKa of said compound in said solution 
having" said anomaly is obtained by determining the pH of 

10 the intersection of horizontal and diagonal asymptotes 
of said solubility-pH curve, 

c. the difference between (1) the pKa of said 
compound in buffer free of anomaly and (2) the pK. of 
said compound in said solution is determined, and 

15 d. said difference is added to the log of the 

apparent intrinsic solubility to determine the log of 
the intrinsic solubility of said compound free of 
anomaly, 

provided that, (D for ionizable acids, pKa of said 
20 compound in said solution minus the pKa of said compound 
in buffer free of anomaly is greater than zero, or (2) 
for ionizable bases, pKa of said compound in said 
solution minus the pK^ of said compound in buffer free 
of anomaly is less than zero, 

14. The method of claim 1 wherein said 
spectrophotometric property is determined at more than 
one point, and the values at said points are corrected 
to account for the presence of interferences, said 
30 corrections being based on 
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a comparing the shape of the spectrophotometric 
property curve of said reference solution to the shape 
of a curve of said compound in saturated- solution, and 

b applying weighting factors, wherein those 
portions of the curve of said reference solution that 
differ from the corresponding ourve prepared 3.n 
saturated solution contribute less to said method than 
those where said shape is similar to said curve relating 
to said saturated solution. 

15. The method of claim 14 wherein said 
spectrophotometric property is OV absorbance. 

16 The method of claim 15 applied to determination of 
15 concentrations of sample solutions when comparing UV 

absorbance spectra of sample solutions to UV absorbance 
spectra of reference solutions. 

17 The method of claim 1, wherein said buffer provides 
20 buffering capacity over the range of pH 3-10 and has 

low-UV absorbance. 

18. The method of claim 17 wherein said buffer 
comprises acetic acid, MES, HEPES and boric acid. 

19 The method of claim 17, wherein said buffer 
con^rises glycolic acid, MES, HEPES, taurine, and 
acetonitrile • 

30 20 An ■ instrument for determining the solubility of a 
compound that is sparingly soluble in water and in a 
buffer, said instrument comprising means for: 



- 71--- 



wo 01/55698 



PCT/USOl/02377 



a. 
b. 



accurately titrating components of a solution, 
collecting components of a solution, 
c mixing solutions, 

d. filtering insoluble residues from solutions, . 

e adjusting pH values of solutions, 

f . determining a spectrophotometric property of a 

solution, and 

performing calculations, includrng these 

relating to said nv absorbance spectra. 

' 21 The mstruaent ot- claim 20 for determining the 
solubility 0£ a compound that is sparingly soluble .n 
„ater and in a buffer, .herein sa.d means for 
determining a spectrophotometric property of a solution 
.5 is selected from the group consisting of a DV range 
spectrophotometer, a visible range spectrophotometer a 
colorimeter, a light scattering detector, a polar.meter 
an optical rotation, a tluori^eter, or a crrcular 
dichroism detector. 

22 A buffer providing buffering capacity over the 
range of pH 3-10 and having low-OV absorbance. 

23. The buffer of claim 22 comprising sodium acetate, 
25 MES monohydrate, HEPES and boric acid. 

24 The buffer of claim 22 comprising glycolic acid, 
MES monohydrate, HEPES, taurine, and acetonitrile. 
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